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ABSTRACT 



o ^ 

fH ^ We present high-resolution spectroscopic observations of GRB 060418, obtained with VLT/UVES. These observations 

were triggered using the VLT Rapid- Response Mode (RRM), which allows for automated observations of transient 
I ' phenomena, without any human intervention. This resulted in the first UVES exposure of GRB 060418 to be started 

O ' only 10 minutes after the initial Swift satellite trigger. A sequence of spectra covering 330-670 nm were acquired at 

^ , 11, 16, 25, 41 and 71 minutes (mid-exposure) after the trigger, with a resolving power of 7 km s~^, and a signal-to- 

C/3 . noise ratio of 10-15. This time-series clearly shows evidence for time variability of allowed transitions involving Fe II 

' fine-structure levels {^1)7/2, ^^5/2, ^03/2, and ^Di/2), and metastable levels of both Fe II (^F9/2 and ''D7/2) and Ni II 

(''F9/2), at the host-galajcy redshift z = 1.490. This is the first report of absorption lines arising from metastable 
levels of Fe II and Ni II along any GRB sightline. We model the observed evolution of the level populations with 
' three different excitation mechanisms: collisions, excitation by infra-red photons, and fluorescence following excitation 

5— ( ' by ultraviolet photons. Our data allow us to reject the coUisional and IR excitation scenarios with high confidence. 

The UV pumping model, in which the GRB afterglow UV photons excite a cloud of atoms with a column density A^, 
distance d, and Doppler broadening parameter b, provides an excellent fit, with best-fit values: log A'^(Fe II) = 14.75lg Q4, 
log N{m II) = 13.84±0.02, d = 1.7 ±0.2 kpc, and 6 = 25±3 km s~\ The success of our UV pumping modeling implies 
that no significant amount of Fe II or Ni II is present at distances smaller than ~1.7 kpc, most likely because it 
is ionized by the GRB X-ray/UV flash. Because neutral hydrogen is more easily ionized than Fe II and Ni II, this 
minimum distance also applies to any H I present. Therefore the majority of very large H I column densities typically 
observed along GRB sightlines may not be located in the immediate environment of the GRB. The UV pumping fit 
also constrains two GRB afterglow parameters: the spectral slope, /3 =— 0.5ti oi and the total rest-frame UV flux that 
irradiated the cloud since the GRB trigger, constraining the magnitude of a possible UV flash. 
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1. Introduction 

The influence of a 7-ray burst (GRB) explosion on its en- 
vironment has been predicted to manifest itself in various 
ways. Strong observational evide nce (.G alama et al. 1998;; 
IStanek et all 120031: iHiorth et af] l2003|) indicates that at 
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* Based on observations collected at the European Southern 
Observatory, Chile; proposal no. 77.D-0661. 



least some GRB progenitors are m assive stars (jWooslevI 
I1993I: iMacFadven fc WooslevI 119990 . and therefore the ex- 
plosion is likely to take place in a star-forming region. 
As the GRB radiation ionizes the atoms in the envi- 
ronment, the neutral hydrogen and metal column den- 
sities in the vicinit y of the explosion are expected to 
evolve with time (Pern a fc Loebl 1X9981 ; IVreeswiik et alj 
I2OOII ; iMirabal et al...2002.) . Ultra-violet (UV) photons will 
not only photo-dissociate and photo-ionize any nearby 
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Table 1. Log of UVES observations 



UT start 


epoch Fig 


. [2] line colour 


AT^ 


'^central 


coverage 


exptime seeing'^ 


FWHM 


S/N 


OTmag'^ 


(2006 April 18) 






(min) 


(nm) 


(nm) 


(min) 


(") 


(km s^^) 


peak 




3:16:07 


1 


black 


11.47 


390 


328-452 


3 


1.5 


6.9 


6 


R (6min) = 14.0 


3:16:02 


1 




11.38 


564 


462-560; 568-665 


3 


1.1 


7.2 


14 


V (6.5min) = 15.0 


3:20:17 


2 


red 


16.61 


390 


328-452 


5 


1.5 


6.9 


6 


z (16.3min) = 14.4 


3:20:12 


2 




16.52 


564 


462-560; 568-665 


5 


1.0 


7.2 


15 




3:26:27 


3 


blue 


25.13 


390 


328-452 


10 


1.5 


6.9 


7 




3:26:17 


3 




24.96 


564 


462-560; 568-665 


10 


1.0 


7.2 


16 




3:38:01 


4 


green 


41.20 


390 


328-452 


20 


1.4 


6.9 


7 




3:37:56 


4 




41.12 


564 


462-560; 568-665 


20 


0.9 


7.2 


17 




3:59:30 


5 


magenta 


70.99 


390 


328-452 


40 


1.4 


6.9 


7 


I (69min) = 16.5 


3:59:24 


5 




70.88 


564 


462-560; 568-665 


40 


0.9 


7.2 


17 


z (78min) = 16.2 


4:41:51 


6 




128.09 


437 


376-498 


80 


1.0 


6.1 


9 


V (lOOmin) = 18.8 


4:41:46 


6 




128.00 


860 


670-852; 866-1043 


80 


0.8 


6.5 


21 


I (135min) = 17.4 



^ Time of flux-weighted mid-exposure since GRB trigger, assuming the light curve decay index a = —1.1. 

^ The seeing has been estimated from the 2-D spectra. 

Approximate magni tud e of the opt ical tr ansient around the time of our spectra, in filters V (jSchadv fc Falcon^ 120061 1. R 
l|Melandri et al.ll2006D . I l|Cobbll2006l ) or z l|Nvsewander et al.ll2006l 1. 



molecular hydrogen, but also quickly excite H2 at larger 
distances to its vibrationally excited metastable levels, 
which can be observed in absorption (jDraine & Hao 2002) . 
Finally, du st grains can be destroyed up to tens of par- 
sees away (IWaxman fc Draind 2000t Fruchter et al. 200X1; 



Draine fc Had l2002t iPerna &: Lazzatil 120021: iPerna et all 
20031) ■ Detection of these time-dependent processes, with 



timescales ranging from seconds to days in the observer's 
frame, would not only provide direct information on the 
physical conditions of the interstellar medium (ISM) sur- 
rounding the GRB, but would also constrain the proper- 
ties of the emitted GRB flux before it is attenuated by 
foreground absorbers in the host galaxy and in interven- 
ing gas clouds. In the X-ray, evidence h as been found 
for a time-variable H I c olumn density (jStarling et al.l 
l2005t ICampana et al.|[2007D . presumably due to the ion- 
ization of the nearby neutral gas. In the optical, none of 
these processes have been obse rved until recently, when 
iDessauges-Zavadskv et al.l ()2006D reported a '-^Sa variabil- 
ity detection of Fe II ^T>7/2 A239(|3, observed at two epochs 
roughly 16 hours apart. Such observations are technically 
very challenging because high-resolution spectroscopy com- 
bined with the rapidly decaying afterglow flux requires im- 
mediate follow-up with 8-10 m class telescopes. 

The Swift satellite, launched in November 2004, has per- 
mitted a revolution in rapid spectroscopic follow-up obser- 
vations, providing accurate (5") positions for the major- 
ity of GRBs within a few minutes of the GRB trigger. 
Numerous robotic imaging telescopes react impressively 
fast (within 10 sec) to Swift triggers. As for spectroscopic 
observations, a number of target-of-opportunity programs 
at most major observational facilities are regularly yielding 
follow-up observations of the GRB afterglow at typically an 
hour after the Swift alert. However, most of these programs 
require significant human coordination between the science 
team and telescope personnel/observers. At the European 
Southern Observatory's (ESO) Very Large Telescope (VLT; 



^ Lines arising from fine-structure levels are sometimes indi- 
cated with stars, e.g. Fe 11* for Fe II '^'Dr/2, Fe 11** for Fe II 
^D5/2, etc.; in this paper we will instead list the transition lower 
energy level term and J value in order to indicate all levels that 
we will discuss in a consistent manner. 



consisting of four unit telescopes of 8.2 m each), a Rapid 
Response Mode (RRM) has been commissioned to provide 
prompt follow-up of transient phenomena, such as GRBs. 
The design of this systerr|3 allows for completely automatic 
VLT observations without any human intervention except 
for the placement of the spectrograph entrance slit on the 
GRB afterglow. The typical time delay, which is mainly 
caused by the telescope preset and object acquisition, is 
5-10 minutes, depending on the GRB location on the sky 
with respect to the telescope pointing position prior to the 
GRB alert. The data presented in this paper are the result 
of the first automatically-triggered RRM activation. 

This paper is organized as follows: the UVES observa- 
tions and data reduction are described in Sect. [21 followed 
by Sect. [31 in which we discuss general properties of the ab- 
sorption systems at the host-galaxy redshift from the detec- 
tion of resonance lines. In Sect.[31 we focus on the detection 
of variability of transitions originating from fine-structure 
levels of Fe II, and metastable levels of Fe II and Ni II. The 
time evolution of the level population of these excited levels 
is modeled in Sect.[5l The results and their implications are 
discussed in Sect. [51 and we conclude in Sect. [71 

2. UVES observations and data reduction 

On April 18 2006 at 3:06:08 UT the Swift Burst 
Alert Telescope (BA T) triggered a 7-ray burst alert 
([Falcone et al.l [2006aD . providing a 3' error circle local- 
ization. Observations with the Swift X-Ray Telescope 
(X RT) resulted in a 5" position about one minute later 
(.Falcone et al.l l2006b[ ) , which triggered our desktop com- 
puter to activate a VLT- RRM request for observations with 
the Ultra-violet and Visual Echelle Spectrograph (UVES). 
This was received by the VLT's unit telescope Kueyen at 
Cerro Paranal at 3:08:12 UT. The on-going service mode 
exposure was ended immediately, and the telescope was 
pointed to the XRT location, all automatically. Several min- 
utes later, the night astronomers Stefano Bagnulo and Stan 
Stefl identified the GRB afterglow, aligned the UVES sht 
on top of it, and started the requested observations at 3:16 
UT (i.e. 10 minutes after the Swift 7-ray detection). This 



see |http://www.eso.org/observing/p2pp/rrm.html| 
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represents the fastest spectral follow-up of any GRB by an 
optical facility (until the RRM VLT/UVES observations 
of GRB 060607, also triggered by our tean i, which were 
start ed at a mere 7.5 minutes after the GRB: lLedoux et all 
|2006() . A series of exposures with increasing integration 
times (3, 5, 10, 20, and 40 minutes, respectively) was per- 
formed with a slit width of 1", yielding spectra covering 
the 330-670 nm wavelength range at a resolving power of 
R = AA/A ^ 43,000, corresponding to 7 km s~^ full width 
at half maximum. These observations were followed by a 80- 
minutes exposure in a different instrument configuration, 
but with the same slit width, extending the wavelength 
coverage to the red up to 950 nm. The dat a were reduced 
with a custom version of the UVES pipeline (jBallester et aL 
I2000D . flux-calibrated using the standard response curved 
and converted to a heliocentric vacuum wavelength scale. 
The log of the observations is shown in Table [1] 

3. Ground-state absorption lines from the host 
galaxy of GRB 060418 

The spectra reveal four strong absorption-line systems at 
redshifts z = 0.603, 0.656, 1.107, and 1.490. In what follows, 
we focus on the highcst-rcdshift absorption-line system at 
2abs = 1-490, which corresponds to the redshift of the GRB 
as shown below; the intervening sy stems are discussed in a 
separate paper (|Ellison et al.ll2006f ). 

At the redshift of the GRB host galaxy, we detect a 
large number of metal-absorption lines, arising from transi- 
tions involving the ground state of various ions (see below), 
fine-structure levels of Si II and Fe II, and metastable lev- 
els of Fe II and Ni II. We will discuss the excited-level lines 
in more detail in Sect. [4l in this section we focus on the 
resonance lines, i.e. lines corresponding to an allowed tran- 
sition from the ion ground state to a higher excited level. 
The ions from which resonance lines are detected are C I, 
C IV, Cr II, Mn II, Si II, Si IV, Zn II, Mg I, Mg II, Ni II, 
Fe II, Al II, Al III, and Ca II. C II A1334 is at the very 
blue edge of our spectrum, where the noise is dominating 
the signal. For Fe I, we determine an upper limit (5cr) on its 
column density of log N(Fe I)<11.48, using Fe I A2484. A 
selection of lines is shown in Fig. [T] Because the resonance 
lines are not found to vary in time, we have combined all 
spectra (see Table [IJ to achieve the highest signal-to- noise 
ratio (S/N) possible. The combined spectrum has S/N=16 
at A=4000 A (Arest=1606 A) and S/N=36 at A=6500 A 
(Arest=2610 A). The redshift corresponding to the zero ve- 
locity has been adopted to be z = 1.49000. At this redshift, 
the Lya line is located at 3027 A, just outside the UVES 
spectral range. 

We now wish to highlight a number of observations that 
can be made from Fig. [1] The vast majority of the column 
density of the low- and high-ionization species (as well as 
fine-structure and metastable species) is located within a 
narrow range of velocity (with a spread of 50-100 km s^^), 
and seems to be contained within two or three main com- 
ponents. This small range in velocit y for the low-ionization 
speci es is also seen in GRB 051111 (iProchaska et al.ll2006l 
|2007|) . Highly saturated lines such as those from C IV, 
Si IV, Mg II, and Al II show components to the blue up to 
—200 km s^^, but these harbour only a small fraction of the 

^ http:/ /www. eso.org/observing/dfo/quality/UVES/qc/ 
response.html 



total column density. The line profile of the high-ionization 
C IV lines follows the profile of the low-ionization lines very 
well, even though the comparison is made difhcult by the 
fact that the C IV lines are m uch stronger. This similar - 
ity is uncommon in QSO-DLAs (jWolfe fc Prochaskal[2000l ). 
The main clump of the low-ionization line profiles, though 
kinematically simple, is a complex mix of broad and nar- 
row components. Thanks to the high signal-to-noise ratio, 
a weak narrow component is clearly observed in the profiles 
of Cr II and Mn II at -20 km s'^. 

For a quantitative analysis, we have simultaneously 
fit Voigt-profiles (using the FITLYMAN context within 
MIDAS) to all resonance lines with at least one non- 
saturated transition. The atomic data required by the 
profile fits (vacuum wavelengths, oscillator strengths and 
damping coefficients) have been taken from Morton (2003i). 
For the oscillato r strengths of N i II AA1709, 1741, and 1751, 
the values from iFedchak et al.i (2000.) have been adopted 
instead. We find that at least three components are re- 
quired to yield an adequate fit to the data. Although a 
3-component fit does not describe the blue side of a few 
high S/N lines perfectly (see e.g. Zn II A2026 in Fig. [J), it 
is the simplest model that fits the data adequately. Adding 
a component on the blue side would also require an ad- 
ditional red component to compensate for the loss of the 
broad component on the red side; the redshift of this addi- 
tional red component would be very hard to constrain. We 
note that the total column density derived would hardly 
change if more components are used. The redshift z and 
the Doppler- broadening parameter h (in km s~^) for each 
component are assumed to be the same for all ions. Only 
for Mg I we had to tweak the redshift of the red component 
for a satisfactory fit. The slightly higher redshift for this 
component is reasonably consistent with the profile of Mg I 
A2852, but it is more probable that the red Mg I component 
is blended. A third possibility is that there are two com- 
ponents in this red feature, where the reddest component 
would have a very high Mg I over Zn II ratio. 

The best-fit redshifts and broadening parameters for 
each component are listed at the top of Table [2l along 
with the fit ionic column densities (individual and the to- 
tal of all components). The fits are shown by the solid 
(red) line in Fig. [T] When comparing the resonance lines 
with the transitions from excited levels of Fe II and Ni II, 
it is apparent that although the profiles are very similar, 
the velocity spread of the latter is smaller. The one ex- 
ception is Al III, whose red component (i.e. the one near 
-1-15 km s~^) is not consistent with the resonance-line fit. 
The colu mn densities that we find are consistent with those 
found by IProchaska et al.l (|2007[ ). with the exception of 
Fe II, where our value of log N(Fe II)=15.07 ± 0.08 is 
lower (at 1.8a significan ce) than their adopted value of 
log N(Fe IIW15.22±0.03. IProchaska et al.l (|2007f) used the 
transitions Fe II A2249 and A2260, which are not saturated. 
In the UVES spectra these lines fall right in the red CCD 
gap, leaving us with one unsaturated line, Fe II A1611, 
which has a lower oscillator strength. Th erefore, we have 
more confidence in the determination by IProchaska et al] 
()2007f ): in the discussion that follows, it should be kept in 
mind that our total Fe II column density is probably too 
low by about 0.15 dex. 

From the column densities we calculate the abundance 
ratios of several ions with respect to iron for each compo- 
nent separately and for the total, adopting the solar val- 
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Table 2. Ionic column densities and abundance ratios in the combined spectrum (see FiglT]for the corresponding profile 
fits). 



Ion Lines used log TV ± aiog jv'' 

component 12 3 total 

2abs 1.49001(1) 1.48999(7) 1.49018(3) 

b (km s"^) 23.6 ±0.4 3.0 ± 0.4 6.5 ±0.1 

CT 1656 12.52 ±0.18 12.29 ± 0.14 12.19 ±0.17 12.83 ±0.11 

Mg I 2026 < 12.07 13.01 ± 0.02 13.60 ± O.Ol'' 13.70 ± 0.02*" 

Sill 1808 15.74 ± 0.02 14.85 ± 0.17 15.34 ± 0.05 15.92 ± 0.03 

CrII 2056,2062, 2066 13.51 ±0.01 12.46 ± 0.06 12.88 ± 0.03 13.63 ± 0.02 

Mnll 2576, 2594, 2606 12.97 ± 0.01 11.82 ± 0.05 12.59 ± 0.02 13.14 ±0.01 

Fell 1611 14.88 ±0.10 14.03 ± 0.31 14.49 ±0.11 15.07 ± 0.08 

Nill 1709, 1741,1751 13.77 ±0.03 12.82 ±0.13 13.03 ± 0.08 13.88 ± 0.03 

Znll 2026,2062 12.80 ± 0.02 11.98 ±0.05 12.69 ± 0.01 13.09 ± 0.01 
Abundance ratio 

[Si/Fe] 0.79 ±0.10 0.75 ± 0.35 0.78 ±0.12 0.78 ±0.09 

[Cr/Fe] 0.45 ±0.10 0.25 ± 0.32 0.21 ±0.11 0.38 ±0.08 

[Mn/Fe] 0.06 ±0.10 -0.24 ± 0.31 0.07 ±0.11 0.04 ±0.08 

[Ni/Fe] 0.14 ±0.10 0.04 ± 0.34 -0.21 ±0.14 0.06 ± 0.09 

[Zn/Fe] 0.76 ±0.10 0.79 ± 0.31 1.04 ±0.11 0.86 ±0.08 



The errors listed are the formal errors provided by FITLYMAN; although an error of 0.02 on an individual component is likely 
to be an underestimate, we consider the error on the total column density to be realistic. 
^ The third component of Mg I is likely to be blended; see the discussion in the text. 



ues from iLoddersI (|2003l ). The resulting values are listed 
in Table ^ The ratio [Zn/F e] is high compared to the 
globa l QSO-DLA population (jLedoux et al.ll2003t IVladild 
I2004D . and suggests a large dust depletion, especially in 
component 3 where [Zn/Fe] = 1.0. The solar value that we 
find for [Mn/Fe] pr ovides additional evidence for substan- 
tial d ust depletion (jLedoux et al.ll20d^ iHerbert-Fort et al] 
|2006() . Given these indications for a large dust depletion, 
the actual value for [Si/Zn] may be 0.2-0.3 dex higher than 
observed ([Si/Zn]tot = —0.08), which would suggest an a- 
element overabundance, provided that zinc can be used as 
a proxy for iron peak elements. Although the values for 
[Zn/Fe] and [Mn/Fe] are high compared to those found 
in QSO-DLAs ([Zn/Fe]Qsos=0-l and [Mn/Fe]Qsos=-0.5- 
0.4), they are rather typical for the ISM of GRB host 
galaxies, with [Zn/Fe]GRBs=l-2 and [Mn/ Fe]GRBs^ 0.1-0 3 
(|Savaglio et al.ll2003l : ISavaelio fc F"aiill2004t Savagho 20061) . 
This dust depletion difference between QSO-DLAs and 
GRB host galaxies can be naturally explained if QSO sight- 
lines on average do not probe the central regions of galax- 
ies (f o r which there is growing evidence, e.g. .Wolfe fc ChenI 
[200a lEllison et al.l [2005t iChen et all l2005ah . while GRB 
lines-of-sight do. 



4. Detection and variability of transitions involving 
excited levels of Fe II and Ni II 

In the left panel of Fig. [2] we show the profiles of some 
selected resonance lines and transitions arising from all 
four fine-structure levels of Fe II (^07/2, ^05/2, ^03/2, and 
^Di/2), as well as from transitions from metastable levels of 
Fe II (4F9/2 and ^07/2) and Ni II (^Fg/z). See Figs.Handli] 
for an illustration of the relevant energetically lower levels 
of Fe II and Ni II, including the first higher excited level, 
and the wavelength and spontaneous decay probability of 
the transitions between the levels. Back to Fig. [21 we over- 
plot the series of five spectra, epoch 1-5 (see Table [TJ, with 
the colours black, red, blue, green and magenta, respec- 



tively. Comparison of the two panels shows clear evidence 
for variability of the excited-level lines, while the strengths 
of the resonance lines are constant in time. To show this 
variability more clearly, we have constructed apparent col- 
umn density profiles based on pixel optical depths in com- 
posite spectra (Savage & Sembach 1991') for each of the 
Fe II fine-structure levels and the Fe II and Ni II metastable 
levels; these profiles, which have been smoothed with a box- 
car of 5 pixels, are shown in the right panel of Fig. [2] 

We estimate the formal significance of this variability 
by measuring the equivalent width (EW) of the individ- 
ual lines in between —30 km s~^ and +40 km s~^ at the 
various epochs, conservatively adding 3% of the EW to its 
formal error, due to the uncertainty in the placement of the 
continuum. Using the different EW values over the differ- 
ent epochs and its mean, we calculate the chi-square and 
the corresponding probability with which a constant equiv- 
alent width can be rejected. For the individual lines shown 
in Fig.H Fe II A2333, Fe II A2607, Fe II A2407, Fe II A2629, 
Fe II A1702, Fe II A2563, and Ni II A2217, the significances 
are 4.5cr, 5.8cr, 2.1cr, 0.3a, 2.5ct, I.Tct, and 3.5ct, respec- 
tively. Using several transitions originating from the same 
level (the same that have been used to construct the appar- 
ent column density profiles in Fig. ^ , we find the following 
numbers: 8.7cr (Fe II ^07/2), 7.4cr (Fe II ^05/2), 3.2cr (Fe II 
6D3/2), 0.5a (Fe II 6D1/2), 2.2a (Fe II 4F9/2), 1.7a (Fe II 
4D7/2), and 2.5a (Ni II 4F9/2). 

We have performed Voigt-profile fitting to the lines orig- 
inating from the excited levels of Fe II and Ni II, inde- 
pendent from the resonance-line fit. The atomic data from 
iMortonI (|200l were adopted when available, and if not we 
have assumed the values from Kurucz (2003) (note that for 
Ni II we have divided the Kurucz oscillator strengths by 
two; see the discussion in Sect. 15. 2p . The Voigt-profile fit 
results are shown in Table [3] and Fig.[3l we only show the fit 
profiles for epoch 3 as the other epochs display very similar 
results, but with somewhat different signal-to-noise ratios. 
Just as with the resonance lines, a satisfactory fit is found 
when using three components. In an initial fit, the redshift 
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Fig. 1. Absorption-line profiles for a variety of transitions detected at the GRB 060418 redshift. To all low-ionization 
species containing at least one non-saturated transition, we have performed simultaneous Voigt-profiles fits using a three- 
component model; the resulting fits are shown by the solid (red) line (see Table[2]for the fit results). The relative velocity 
of the different components are indicated by the (blue) vertical dotted lines; note that two components have very similar 
redshifts. It is clear that this 3-component fit does not describe some high S/N lines, such as Zn II A2026, very well; 
we will come back to this in the discussion. We note, however, that the total column density would hardly change if 
additional components would be introduced. For comparison purposes, we also show transitions that we have not fit: 
(saturated) higher-ionization lines of C IV and Si IV, lines from Fe II fine-structure levels, and transitions originating 
from metastable levels of both Fe II and Ni II. 



and b parameter were free to vary from epoch to epoch. As 
these were found to be constant with time, in a final fit they 
were fixed to the averages over the five epochs; the z and 
6 averages and the corresponding standard deviations are 
listed at the top of Tabled The column density errors listed 
in Table [3] are the formal errors provided by FITLYMAN. 



We have also estimated an error in the placement of the 
continuum by varying the sigma clipping factors and the 
order of the polynomial with which we fit the continuum 
around each line, and rerun the Voigt profile fit. The max- 
imum change that we find is 0.03 dex, which we add to the 
formal error for the rest of the analysis. 
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Fig. 2. The epoch 1-5 UVES spectra of GRB 060418 (see Tabled]) are overplotted with the colours black, red, blue, green 
and magenta, respectively. In the left panel individual lines are shown, typical resonance lines on the left, and the lines 
arising from the excited levels of Fe II and Ni II on the right. The latter show evidence for a varying equivalent width as 
a function of time. To make this variability clearer, we have combined various lines that arise from the same level and 
constructed apparent column density profiles, smoothed with a boxcar of 5 pixels; these are shown in the right panel. 



iProchaska et al.l (j2007[ ) have performed time-resolved 
high-resolution spectroscopy of GRB 060418 as well. Their 
three spectra were taken around the same mid-exposure 
time as our epoch 4, 5 and 6 spectra. They do not con- 
sider variability, and report on the average column density 
for the four Fe II fine-structure levels; they do not mention 
the metastable levels of Fe II and Ni II. When comparing 
their average values with our epoch 4 column densities, the 
results are fully consistent within the errors. 



Comparison of the resonance-lines fit with the excited- 
lines fit shows that the redshifts and b parameters for the 
three components are very similar, cf. Tables [5] and [31 When 
we run the excited-lines profile fit with the redshift and b 
parameter fixed at the values of the resonance lines, the re- 
sulting fit is very poor. Therefore, although the fits provide 
similar results, the redshifts and b parameters are signifi- 
cantly different. We already noted this difference in Sect. [31 
and we will come back to this point in Sect. [51 
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Fig. 3. Absorption-line profile fits to selected transitions from excited Fe II and Ni II. The lower level of the transition, 
for which the column density is determined from the fits, is indicated in each panel. The fit results are listed in Table [31 
We only show the fits for the third epoch spectra; the quality of the spectral fits for the other epochs are very similar, 
with only small differences due to slightly different signal-to-noise ratios (see Table [T]). 



Atomic fine-structure levels are caused by an energy 
split due to the interaction of the total electron spin and 
total angular momentum of the electrons. The transitions 
between these levels are not allowed, i.e. they cannot pro- 
ceed through an electric dipole transition, and therefore the 
corresponding transition probabilities are low. The same is 
applicable to other energetically lower levels, also called 
metastable levels, and their fine-structure levels. Figures |3] 
and [5] show the energy level diagrams, of selected levels of 
Fe II and Ni II. 

These fine-structure and metastable levels can be pop- 
ulated through (1) collisions between the ion and other 
particles such as free electrons, (2) direct photo-excitation 
by infra-red (IR) photons (with specific wavelengths be- 
tween 87-260 /im), and/or (3) indirectly through excita- 
tion by ultra-violet (UV) photons, followed by fluorescence. 
Detection of transitions from these energetically lower ex- 
cited levels provides a powerful probe of the physical con- 
ditions in the interstellar medium (Bahcall & Wolf 1968), 
where the quantities that can be derived depend on the 
excitation mechanism. 



IVreeswiik et al.l (j2004l) noted the presence of transi- 
tions originating from the fine-structure level of Si II in 
the host galaxy of GRB 030323; as these lines had never 
been clearly detected in QSO-DLAs (we note that they 
had been dete cted in absorption systems associated with 
the Q SO, see IWampler et al.l 119951 : ISrianand fc PetitjeanI 
|2001[ ). this detection suggested an origin in the vicinity 
of the GRB. Assuming that collisions with electrons were 
the dominant excitation mechanism, a volume density of 
uhi = 10^ — 10^ cni^'^ was derived (see also Savaglio & Fall 
I2OOI IFvnbo et al.H200l . We note that Si II 2P3/2 A1816 
is not detected (5 a upper limit: log N< 14.81) in the case 
of GRB 060418, and that Si II ^Pg/a A1533 (see Fig. H]) is 
severely blended with Fe II A2382 at the redshift of an in- 
tervening absorber, z — 0.603. Therefore, this (or any) Si II 
fine-structure level is not included in our analysis. 

More recently, even more exotic transitions involv- 
ing fine-structur e levels of Fe II have been discovered in 
GRB sightlines JChen et all l2005bt iPenprase et"!!! 120061: 



Prochaska et al 



Prochaska et al 



20061: ID 'Elia et al.ri2006D . As noted by 
( 20061 ). these lines had previously been 
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Table 3. Column densities of the Fe II and Ni II excited 
levels, in the three individual components and the total, at 
the various epochs (see Table [IJ since the burst trigger. 



Fe II 



epoch 




log Af±aiog jv"" 


comp. 


1 


2 3 total 




1.48996(6) 


1.49002(4) 1.49015(9) 


b 


28.0 ±3.5 


5.0 ±2.1 10.2 ±1.5 

6T^ NNnooo nnon nnon n^-io r,^r,^ 



13.49 ±0.02 
13.52 ±0.02 
13.48 ± 0.02 
13.37 ±0.02 
13.23 ±0.03 



J7/2 ■ 
13.27 ±0.03 
13.17 ±0.03 
13.14 ±0.03 
13.13 ±0.03 
13.06 ±0.03 



13.63 ± 0.02 

13.64 ±0.02 
13.62 ±0.02 
13.57 ±0.02 
13.41 ±0.02 



13.96 ±0.01 
13.96 ±0.01 
13.93 ±0.01 
13.87 ±0.01 
13.73 ±0.01 



Fe irD5/2 AA2328, 2381, 2399, 2607, 2618 



13.09 ±0.05 
13.12 ±0.04 
13.12 ±0.04 
13.04 ± 0.04 
12.84 ±0.07 



J5/2 
13.04 ±0.03 
13.06 ±0.03 
13.02 ±0.03 
12.96 ±0.03 
12.86 ±0.04 



13.45 ± 0.02 
13.45 ± 0.02 
13.39 ± 0.02 
13.36 ±0.02 
13.24 ±0.02 



13.71 ±0.02 

13.72 ± 0.02 
13.68 ±0.02 
13.63 ±0.02 
13.50 ± 0.02 



Fe IPDg/a AA2338, 2359, 2407, 2411, 2614, 2621 



^3 /2 

12.60 ± 0.09 
12.63 ± 0.08 
12.59 ±0.08 
12.56 ±0.08 
12.38 ±0.13 



12.88 ±0.03 

12.85 ±0.03 

12.86 ±0.03 
12.78 ±0.03 
12.67 ±0.04 



13.17 ±0.02 

13.18 ± 0.02 
13.18 ±0.02 
13.13 ± 0.02 
13.00 ±0.02 



13.42 ±0.02 
13.42 ± 0.02 
13.42 ± 0.02 
13.36 ±0.02 
13.23 ±0.03 



Fe IPDi/a AA2345, 2411, 2622, 2629 



12.38 ±0.15 

12.39 ±0.14 
12.65 ±0.07 
12.44 ±0.10 
12.35 ±0.13 



1/2 

12.51 ±0.06 
12.42 ±0.06 
12.34 ±0.07 
12.42 ±0.05 
12.21 ±0.08 



12.77 ±0.03 
12.80 ±0.03 
12.77 ±0.03 
12.79 ±0.03 
12.73 ±0.03 



13.06 ± 0.04 
13.06 ± 0.04 
13.10 ±0.03 
13.06 ± 0.03 
12.97 ±0.04 



Fe II^Fg/a AA1566^ 1612^ 1637^ 1702^ 2332, 2360 

1 12.74 ±0.33 12.79 ±0.14 13.17 ± 0.06 13.42 ±0.10 

2 12.99 ±0.17 12.62 ±0.19 13.26 ± 0.05 13.51 ±0.07 

3 12.91 ±0.19 12.89 ±0.10 13.40 ± 0.04 13.61 ±0.05 

4 12.97 ±0.15 13.09 ±0.06 13.41 ± 0.03 13.68 ± 0.04 

5 12.54 ±0.44 13.03 ±0.07 13.47 ± 0.03 13.64 ± 0.06 

6 13.38 ±0.12 12.53 ±0.42 13.42 ± 0.07 13.73 ± 0.08 
FT? 



Fe IPD7/2 AA1635, 2563 



12.36 
12.33 
12.46 
11.35 
12.45 



±0.30 
±0.30 
±0.20 
±0.80 
±0.20 



12.62 
12.83 
12.85 
12.82 
12.86 



Ni II 
±0.11 
±0.07 
±0.06 
±0.06 
±0.06 



J7/2 
12.35 ±0.15 
12.09 ±0.25 
12.35 ±0.12 
12.41 ±0.10 
11.54 ±0.72 



12.78 ±0.06 
12.75 ± 0.06 
12.40 ± 0.12 
12.49 ± 0.09 
11.98 ±0.30 



13.02 
12.95 
12.88 
12.77 
12.61 



±0.10 
±0.11 
±0.10 
±0.10 
±0.20 



'9/2 



AA2166, 2217, 2223, 



12.59 ±0.06 
12.63 ±0.06 
12.66 ±0.05 
12.68 ±0.05 
12.61 ±0.06 



13.03 ± 0.03 
13.10 ±0.02 
13.12 ±0.02 
13.22 ±0.02 
13.24 ±0.02 



2316 
13.27 
13.37 
13.40 
13.45 
13.46 



±0.03 
±0.03 
±0.02 
±0.02 
±0.02 



The errors listed are the formal errors provided by 
FITLYMAN; although an error of 0.02 on an individual 
component is likely to be an underestimate, we consider the 
error on the total column density to be realistic. 

These lines are also covered in the spectrum with setting 
437 nm (see Table [T|, resulting in the determination of column 
densities at a 6* epoch. 



detected in absorption in extreme envi ronments such a s 
Broad Abs orption-Line (BA L) quasars (|Hall et al.|[200l . 
rj Carinae (iGuU et al.l I2005D . and the disk of /3 Pictoris 
(|Lagrange-Henri et al.l 1198 8*). For the Fe II fine-structure 
level population along GRB sightlines it has been argued 
(jProchaska et al.| |2006[ ) that IR excitation is negligible, 
that collisional excitation is improbable (although not ex- 
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Fig. 4. Energy level diagram for selected levels of Fe II. For 
the lower levels we only show the levels for which we detect 
transitions, i.e. the fine-structure levels of the Fe II ground 
state, and ^Fg/2 and ^07/2- Note that for clarity reasons, 
we do not show the fine-structure levels of the latter. With 
the arrows, we indicate the most likely transitions between 
these levels, including one higher excited level. For each 
transition we show the wavelength in A and the sponta- 
neous decay Einstein coefficient, Aui in (which is pro- 
portional to the absorption coefficient Biu, see also below). 
The electric dipole allowed transitions are indicated with a 
solid line, and the forbidden transitions (magnetic dipole or 
electric quadrupole) with a dotted line. Note that to pop- 
ulate the level either four IR photons are required, 
or two UV photons, where the higher levels involved need 
to have J=7/2 and 3=3/2. This, combined with the much 
larger transition probabilities in the UV, makes the UV 
pumping mechanism much more efficient than IR excita- 
tion. 



eluded), and that indirect UV pumping probably is the 
dominant excitation mechanism. The detection of variabil- 
ity at the 3cr level (using two different instruments) of one 
Fe II fine-stru cture line in the spectrum of G RB 020813 
was reported (|Dessauges-Zavadskv et al.ll2006[ ). which the 
authors claim to be supportive evidence for the UV pump- 
ing model. 

The lines arising from the metastable levels of both 
Fe II C^Fg/^ and ''Dr/a) and Ni II C^Fg/a) that we detect 
are the first lines from metastable levels to be identified 
along any GRB sightline. However, t hese have also been 
previously detected in BAL quasars (iHazard et al. 1987t 
IWampler et"allll995f ) and rj Carinae (|Gull et alJl2005l )~ 
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Fig. 5. Energy level diagram for selected levels of Ni II (see 
Fig. [4]). Note the very low A- value for the transition '^Fg/2 
to the ground state (this is actually an electric quadrupole 
transition), while the level ^Fg/2 can be very easily pop- 
ulated from the higher excited levels, one of them being 
^Dy^2- Therefore this level can be expected to be densely 
populated in the presence of a strong UV radiation field. 
The fine-structure level of Ni II ground state, ^03/2 has 
a relatively high probability of spontaneous decay, with a 
mean lifetime of 1/(5.4x10^^ s"^)=19 s. 



Our clear detection of time- variation of numerous tran- 
sitions involving all fine-structure levels of the Fe II ground 
state, and moreover transitions originating from metastable 
levels of Fe II and Ni II, allows for a critical comparison 
of the data with the three possible excitation mechanisms 
mentioned above. However, independent of the mechanism 
at play, the detection of time-variable absorption implies 
that the flux from the GRB prompt emission and/or after- 
glow, directly or indirectly, is the cause of the line variabil- 
ity, and that the absorbing atoms are located in the relative 
vicinity of the GRB explosion. 



5. Modeling of the time evolution of Fe II and 
Ni II excited levels 

5.1. Collisional excitation 

We first cons i der t he collisional model. Although 
iProchaska et al.l (|2006f ) did not detect any change in col- 
umn density of excited Fe II, they did suggest that vari- 
ability of absorption-line strengths would be inconsistent 
with a collisional origin of the excitation. This is certainly 
expected to be the case for a medium out of reach of the 
influence of the GRB afterglow fiux, but close to the GRB 
site one might expect the incidence of intense X-ray and 
UV radiation to deposit a considerable amount of energy in 
the surrounding medium through photo-ionization, causing 
a situation similar to photo-dissociation regions (PDRs), 
which can produce a shock front with typical velocities 
of 10-20 km s"'^ and density enhancements. Relaxation of 



these high-density regions might then result in a change 
in column density. The profiles shown in Fig. [T] are actu- 
ally suggestive of such a shock front situation: they show 
two main components with a velocity difference of 25 km 
s~^, and moreover, the lines caused by transitions from the 
metastable levels seem to be enclosed by the ground state 
species of e.g. Cr II. Therefore, it is quite reasonable to 
consider the collisional excitation scenario. 

If collisions of the Fe II ions with electrons, protons 
or HI atoms is the dominant excitation process, and the 
collisional de-excitation rate exceeds the spontaneous decay 
rate, the population ratio between two lev els i and 7 should 
follow the Boltzmann distribution (see e.g. IProchaska et al.l 
I2OO6D : 



(1) 



where g is the statistical weight of the level, Eij is the 
energy jump between the levels, k is the Boltzmann con- 
stant, and Tex is the excitation temperature. We fit this 
Boltzmann distribution to all available excited levels of 
Fe II at each epoch separately. For each epoch we fit two 
parameters: the density (in our case column density) and 
the excitation temperature. The resulting Boltzmann fit is 
shown in Fig. [51 Even with 2 free parameters for each of 
the five epochs, the model is clearly not able to reproduce 
the observed column densities: the reduced chi-square is 
xl = 95.7/(5-2) = 31.9. 

If we would have only detected the fine-structure lev- 
els of Fe II, the collisio nal excitation model fit would have 
been acceptable, as in IProchaska erall (|2006t) . The mam 
culprit for the poor fit is the increasing level population of 
the metastable level Fe II ''Fg/2, which cannot be accom- 
modated in the Boltzmann fit because all the other levels, 
with similar energies, are decreasing with time. Inclusion 
of Ni II ^Fg/a (E=8393 cm"!) would make the fit even 
worse, as its energy level is similar to the Fe II "^07/2 
level (E=7955 cm~^), while its observed column density 
is increasing with time (sec the bottom panel of Fig. 
Moreover, the best-fit Boltzmann model predicts column 
densities for the fine-structure levels of Fe II ^Fg/2 (e.g. the 
predicted column densities for its first fine-structure level, 
Fe II ^F7/2, is log N — 13.4 — 13.6) that are inconsistent 
with the upper limits we obtain for this level (see Fig. [9]) . 
To preserve clarity, we do not show these in Fig. jS) 

An implicit assumption in this collisional excitation 
model is that of local thermodynamic equilibrium (LTE), 
while the observed variability of the absorption lines sug- 
gests that this may not be valid. However, using the 
PopRaticH code (jSilva fc Viegasll200l . we find that if col- 
lisions is the dominant excitation mechanism, the observed 
population ratios of the Fe II fine-structure levels require 
an electron volume density of at least rtg ~ 10'* cm""^. As 
this is very high, while at the same time the observed vari- 
ability is relatively smooth in time, the assumption of LTE 
probably is valid. 

In conclusion, the collisional model is rejected with high 
confidence. 



* We note that PopRatio only includes collisions with elec- 
trons, which are the dominant collision partners for tempera- 
tures below approximately 100,000 K; beyond this temperature 
the contribution of collisions wi th protons and HI ato ms become 
significant (see Figs. 3 and 7 of lSilva fc Viega3l2002l ) 
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Fig. 6. The top panel shows the observed total column den- 
sities (see the last column of Table ^ for the fine-structure 
lines (open circles; from top to bottom: ^T)7/2, ^Ds 721^03/2, 
and ®Di/2, respectively), the first metastable level (filled 
triangles, ^Fg/2), and the second metastable level (filled 
squares, ''D7/2) of Fe II. Overplotted are the results of 
the best- fit Boltzmann model (collisional scenario): solid 
lines for the fine-structure levels, dashed line for '*Fg/2, and 
dashed-dotted for ^Di/2- The best- fit Fe II ground state 
column density is shown by the dotted line, while the best- 
fit excitation temperature (Tgx) is depicted in the bottom 
panel. It is clear that the Boltzmann model can be rejected 
with high confidence. 



5.2. Radiative excitation by GRB-afterglow pliotons 

To verify if our observations can be explained by radiative 
excitation (by IR and/or UV photons), we now consider a 
model of a cloud with column density N (atoms cm~^), at a 
distance d (pc) from the GRB. The afterglow flux will excite 
the atoms in the cloud, and we will calculate the atom level 
populations as a function of time, to be compared with our 
observations. We will only consider excitation, and neglect 
ionization, which, as we will see below, is fully justified. 

We can describe the afterglow flux in the host-galaxy 
rest frame by: 



7rcst 



1.192x10-25 





~/3 


1.083x10^" pc' 


2 


.5439 A. 




d 





in erg s^^ cm^^ Hz^^, where we have used the V=14.99 
UVOT m easurement at 393 seconds after the burst 
(jSchadv fc Falcone 20Qg), corrected for foreground extinc - 
tion in the Galaxy with Ay^O.lA (jSchlegel et all 119981 ). 
and in the absorber at z ~ 1.1 (which shows a clear 
2175 A extinction bump) with a Milky Way extinction 
curve and Ay = 0.25 (at z = 1.1), r esulting in an effec- 
tive Av,z=o — 0.55 (Ellison et al.ll2006l ). So the constant in 
Eq.[2]is the would-be observed UVOT flux at z = if there 
would not have been any foreground extinction. The best- 
fit afterglow intrinsic spectr al slope assum ing these extinc- 
tion values is /3 = -0.8 (see lEllison et al.ll2006f) . However, 
as this value is quite uncertain we will also determine a 
best- fit value for /3 in the fit. The flux decay in time of 
our spectra is very well described by a power law with in- 
dex — 1.1, which we adopt for a. The value for the decay 
index determined from UBVRIz p hotometry data reported 
in GCNs range from -1.1 to -1.3 (iNvsewander et al.|[2006l : 
iCobbll2006l : ISchadv fc Falconell2006f ). For the calculation of 
the luminosity distance to the GRB, d/ = 1.083 x 10^° pc, 
we have adopted Hq — 70 km s"^ Mpc^^, flu = 0.3, and 
r^A = 0.7. 

The atom level population of an upper level u with re- 
spect to a lower level I is given by the balance equation: 



dN, 



dt 



(3) 



where Aui, Bui and Bju are the Einstein coefficients for 
spontaneous decay, stimulated emission and absorption, re- 
spectively, with Bui = Au\X^/2hc (all in cgs units), and Biu 
= Bui gu / gi (g is the statistical weight of the energy level, 
with g=2J-|-l, and where J is the total angular momentum 
of the electrons). Fi/(to) is the incoming afterglow flux at 
the monochromatic frequency corresponding to the transi- 
tion energy, an d modified by t he optical depth at line center 
(see Eq. 3.8 of iLequeuxl [20051) : 



(4) 



with To = i-49i-xiu — ™id and where the oscillator strength 
/ is calculated from Aui, using: 



/ 



(5) 



(2 



Si, is the source function (see Eq. 3.6 of lLeaueiu3 12005[ ) . 
The Doppler width, or broadening parameter, 6, has been 
determined from the line-profile fits to be 28 km s"^, 5 km 
s"^, and 10 km s^^ for the three different components (see 
Table [3]). We allow the b value to vary with the aim to 
obtain the best-fit value, to be compared with the above 
measurements. For many UV transitions the cloud that we 
model will be optically thick, and therefore we slice up the 
cloud in a sufficient number of plane-parallel layers, so that 
each layer can be considered optically thin for a particular 
transition; we set the maximum allowed optical depth for 
a layer to be riayormax=0.05. 

An important ingredient in the model fit is the adopted 
atomic data values for the spontaneous decay coefficients 
Aui (or equivalently, /, see Eq. [5]), and that these are ex- 
actly the same as used in the Voigt profile fits performed 
to obtain the observed column densities (see Tables [5] and 
|3]). We made sure that this is the case. For Fe II we in- 
) elude the 20 lower energy levels in our calculations (up to 



Vreeswijk, Ledoux, Smette et al.: Rapid-Response Mode VLT/UVES spectroscopy of GRB 060418 



11 



E=18886.78 cm~^), a nd the A's b e tween all these lower 
levels are taken from iQuinet et al.l ()1996[ ). For the tran- 
sitions between the lower and hi gher excited levels, we 
adopt the values bv iMortonI (l2003l) if availab le, and if not 
then we use those provided bv lKuruc3 (|2003( !FI. The num- 
ber of Fe II higher excited levels included is 456, with a 
resulting number of transitions of 4443. For Ni II we in- 
clude the lo wer 17 energy levels, and take the A's between 
these from iQuinet fc Le Dourneuj ( 1996D. complemented 
by those from iNussbaumer fc StorevI (|l982l ). For the Ni II 
ground state transitions corre sponding to Ni II A131 7 and 
1370 we adopt the / values of lJenkins fc Tr ipp (200 6|), and 
for Ni II AA1454, 1709, 1741 and 1751 from Fcd chak et all 
(|2000t ). For the other transition probabilities between the 
lower and higher excited levels of Ni II we again use the 
value from [M orton (2003) if available, and otherwise those 
from iKurucj p003f ). As the ratio of the /-values of the 
Jenkins et al. and Fedchak et al. ground-state lines com- 
pared to the Kurucz values varies from 1.87 to 2.59, and 
we find similar ratios between values of two Morton Ni II 
fine-structure lines and those of Kurucz, we have divided all 
Kurucz Ni II A's by a factor of two. We stress that although 
this factor of two results in different inferred column den- 
sities, it does not affect the fit results since we use exactly 
the same oscillator strengths in our model. The number of 
Ni II higher excited levels included is 334, with a resulting 
number of transitions of 3136. 

We have written an IDL routine that incorporates the 
equations above and the adopted Fe II and Ni II atomic 
data values, and calculates the level evolution of the atoms 
in the cloud as a function of time. This model is fit to the 
observations (using Craig Markwardt's MPFIT routine^) 
with the following free parameters: the distance d, the total 
Fe II or Ni II column density N, the afterglow spectral 
slope /?, the Doppler parameter b, and the rest-frame time 
at which we start the calculations, to- We note that this to 
does not provide any constraints on the shape of the light 
curve before the time that our first spectrum was taken 
(we simply extrapolate the light curve back to to assuming 
a decay index of a = —1.1), but it does constrain the total 
number of photons that arrived at the cloud since the GRB 
trigger. 

By selecting the levels that we loop through, we can 
either treat IR excitation and UV pumping separately, or 
combine the two in a consistent manner. For the UV tran- 
sitions, we assume that the higher excited levels are merely 
a route to any lower level that the higher level can com- 
bine with, i.e. after excitation of a number of atoms in 
one timestep, all electrons are re-distributed immediately 
among all possible lower levels and no electrons stay in 
the higher excited level. This is justified by the very large 
spontaneous decay transition probabilities for nearly all 
higher excited energy levels. As a consistency check, we 
compa red the results of o ur program with the PopRatio 
code (Silva fc Viegasl[2002f ). which computes the Fe II fine- 
structure level population assuming an equilibrium situa- 
tion, i.e. = in Eq.[3l Using exactly the same Galactic 
UV background as in PopRatio (converted to flux density 
by applying a factor of 4 tt), and in the optically thin limit, 
the results are identical down to the 0.07% level. 



see 
see 



http: / /kurucz.harvard.edu 
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Fig. 7. The top panel shows the same as the top panel of 
Fig. El but now with the IR excitation model overplotted: 
solid lines for the fine-structure levels, dashed line for '*Fg/2. 
The 12 level fit column density does not even reach the 
lower limit of the plotting range. The model prediction for 
the evolution of the Fe II ground state column density is 
shown by the dotted line. 



5.2.1. IR excitation only 

First we consider only IR photons to be exciting the atoms 
in the cloud. The 20 lower energy levels of Fe II are in- 
cluded, and we do not consider Ni II. The resulting fit is 
shown in Fig. [T] We note that we have imposed a lower 
limit to the distance of 2 pc and we fixed the spectral 
slope ai (3 — —0.8 and the h parameter at 18 km s~^; 
the value for the latter is unimportant as all IR transi- 
tions are basically optically thin. For distances lower than 
2 pc, the calculation would take too long to compute on 
our workstation. The reason for this is that we adjust the 
program timestep in such a way that a maximum of 5% 
of all atoms can be excited to the higher excited level of 
a particular transition in each timestep; for a large photon 
flux this requires a very small timestep, i.e. a large number 
of calculations. If we would allow the distance to go un- 
der 2 pc, the fitting routine would try to move the levels 
^Fg/2 (dashed line) and ^07/2 (below the lower limit of the 
plotting range, peaking at log N = 11.6) up, which would 
also cause the fine-structure levels (solid lines) to move up 
slightly. The final chi-square would be lower than for the 
present > 2 pc fit, which has (IR) = 2571/(31-3) = 91.8, 
but it would still provide an extremely poor fit to the ob- 
servations. Moreover, at such a short distance, most of the 
Fe I I would be expected to be ionized in the first place 
(e.g. lWaxman fc Draind 120001 : iPerna fc Lazzatill2002D . and 
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therefore we can safely reject IR excitation mechanism as 
the explanation for the observed level population and evo- 
lution. 

The reason for the relatively low population of the 
metastable levels compared to the Fe II ground state fine- 
structure levels in the IR excitation case is not due to a 
lower transition probability for the former; e.g. between 
the ground state and its first fine-structure level ^07/2, 
A=2.13xl0~'^ s~^, while between the ground state and 
the second metastable level '^07/2, A=4. 74x10"'^ s~^ (see 
Fig. 2]). The reason is the wavelength dependence to the 
third power of the Einstein absorption coefficient B^, (see 
below Eq. [3]): photons with a longer wavelength are much 
more likely to be absorbed. For the levels mentioned above, 
this makes the transition from the ground state to ^07/2 
a factor of (7955/385)3 ^ gggo less likely, while the differ- 
ence in the observed column density is only a factor of 10. 
Had we only observed the variation of the fine-structure 
levels of the ground state, and not the levels ''Fg/2 and 
'*D7/2, we would have not been able to reject the IR ex- 
citation model with such high confidence, as merely con- 
sidering those levels results in an excellen t fit to the data, 
with Y ^/IRflievciJ = 11.0/(20 - 3) = 0.65. iProchaska erall 
(|2006f ) rejected the IR excitation scenario on the basis that 
IR pumping is negligible at the distance limit set by the de- 
tection of Mg I in their spectra (which assumes that Mg I 
and the excited material is at the same location, which 
need not be the case; see also Sect. combined with the 
observation that UV pumping is dominant at any given dis- 
tance from the GRB, in the absence of severe extinction. 
Although these arguments are strong, they are not as con- 
clusive as our modeling results. 



5.2.2. UV pumping 

After rejection of coUisional and IR excitation, we now con- 
sider the UV pumping scenario. In the UV model calcula- 
tions we consider 20 lower and 456 higher excited levels of 
Fe II. The resulting fit is shown in the top panel of Fig. [S] 
The best-fit values for the fit parameters are as follows: 
log 7V(Fe II ground state)=14.75l^;2t d ^ 1.7 ± 0.2 kpc, 
P =-0.5t?;^, to=74li^ s, and 6 = 25 ± 3 km s-^, and 
a chi-square of x^(UV-Fe II) 26.2/(31 - 5) = 1.01. 



^rest (minutes since GRB trigger) 



Next, we also model the evolution of the Ni II ''Fg/2 level, 
using 17 lower and 334 higher levels of Ni II. We fix all 
parameters in the Ni II fit to the best-fit values of the 
Fe II fit, except for the Ni II ground state column density. 
The resulting fit is shown in the bottom panel of Fig. [51 
The reduced chi-square is xUW ~ Ni II) = 5.6/(5 - 1) = 
1.4, and the best-fit Ni II column density is log A^(Ni II 
ground state)=13.84 ± 0.02. When also including the dis- 
tance as a free parameter, we find log iV(Ni II ground 
state)=13.73 ± 0.02, and d = 1.0 ± 0.3 kpc, with a chi- 
square of x^(UV - Ni II) = 0.72/(5 - 2) = 0.24. 

From Figs, m and [51 it is straightforward to see why the 
levels Fe II ''Fg /2 and Ni II ''Fg /2 increase with time in the 
UV pumping scenario. The route to these levels is rather 
quick: one out of every 5000 photons at 2600 A will bring 
the ion to this excited level. We note that the higher excited 
level shown in Fig. [His just one out of many levels that al- 
low population of the Fe II ^Fg/2 level through absorption 
of a UV photon, followed by spontaneous decay. Once in 
this level, it takes 1/(9.2x10"^ s~-^) = 3.0 hours for the ion 
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Fig. 8. The top panel shows the same as the top panel 
of Fig. [SI but now with the UV pumping model overplot- 
ted: solid lines for the fine-structure levels, dashed line for 
'*Fg/2, and dashed-dotted for '*D7/2- The bottom panel dis- 
plays the observed total column densities for Ni II '*Fg/2 
(filled circles), and the best-fit Ni II model. In this Ni II 
fit, all parameters except for Ni II column density were 
fixed to the best-fit values obtained from the Fe II fit. The 
model prediction for the evolution of the Ni II ground state 
column density is shown by the dotted line. All Fe II and 
Ni II column densities are very well described by the UV 
pumping model. 



to decay to the Fe II ground state; this is longer than the 
time scale over which our spectra were recorded (1 hour in 
the rest frame) , and explains why this level continues to rise 
in Fig. [51 Ni II "^Fg /2 is even easier to populate through the 
absorption of UV photons, and will take a longer time to de- 
cay to the Ni II ground state: 37 hours. Transitions arising 
from these Fe II and Ni II metastable levels are therefore 
excellent probes of the UV pumping mechanism, as they 
can be observed up to many hours after the GRB event. 
In fact, although we were the first to identify them, these 
lines should also be pres ent in the high-resolut i on sp ectra 
of GRBs 050730, 051 111 (IProchaska et al.l l2006l I2007D and 
050922C (|D'Elia et al. l l2005fl . 

As we calculate the level population of the lower 20 
levels for Fe II, and lower 17 levels for Ni II, we can compare 
if the model predictions for all levels are consistent with our 
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Fig. 9. Comparison of the UV-model predicted column 
densities for several additional Fe II and Ni II excited levels, 
i.e. others than the ones shown in Fig. [51 with the Scr upper 
limits that we obtain from the spectra. The Fe II "^1)5/2 level 
is actually detected at just above 5a in the second epoch. 



data. Searching for the detections of lines originating from 
these levels has resulted in one new detection (epoch 2 for 
level Fe II but the rest we can only place an upper 

limit (we adopt 5a) to the column density, as shown in 
Fig. [3 

iDessauges-Zavadskv et al.l (|2006D have reported a sig- 
nificant (~3(t) decline by at least a factor of five in the 
equivalent width of the Fe II ^07/2 A2396 transition, in 
spectra taken at 4.7 and 20.8 hours after GRB 020813. We 
note that this line is saturated in our spectra, and moreover 
blended with Fe II ^05/2 A2396, and therefore we do not 
use it in our analysis. To verify if this decline is roughly 
consistent with our calculations for GRB 060418, we de- 
termine the best-fit UV pumping model Fe II ^07/2 col- 
umn densities for GRB 060418 at 2.1 and 9.2 hours in the 
rest frame, using t he redshift of GRB 020813, z = 1.255 
jBartFetaDIlM). We obtain log N(trcst=2.1 hr)=13.16 
and log N(trest=9.2 hr)=12.45, corresponding to a de- 
cr ease of a factor of 5.1, fu ll y con sistent with the result 
of iDessauges-Zavadskv et aD (|2006( ). 



6. Discussion 

Comparison of our UVES data with modeling of coUisional 
and radiative excitation clearly shows that UV pumping 
by GRB afterglow photons is the mechanism responsible 
for the population of the Fe II and Ni II excited levels. 
However, in the UV pumping model, one would expect the 
ions in the ground state and the excited levels to be ex- 
actly at the same location and velocity. As we discussed 
in Sects. [3] and [H this is not exactly the case. One sim- 
ple way of resolving this apparent discrepancy is to in- 
voke more components than the three that we resolve. A 
fourth and fifth unresolved component could be hidden in 
the resonance-line profile at the same velocities as the two 
main components of the excited levels. Then the observed 
velocity offset between the excited levels and the resonance 
lines can be explained if the observed resonance-line compo- 
nents are mainly due to gas that is further away from the 
GRB, so excitation of the ground state by UV pumping 
is negligible. This requirement of additional ground-state 
components is consistent with the rather poor fit of our 3- 
component Voigt-profile fit to the blue wing of a few high 
S/N resonance lines (see Fig. [T]). The existence of addi- 
tional components is fully consistent with the UV pumping 
fit results. According to this fit, the Fe II and Ni II ground 
state column densities are 0.3-0.5 dex lower than what we 
measure in the spectra, providing evidence for gas along 
the line of sight in the host galaxy that is not affected by 
the GRB, and thus this gas needs to be further away than 
the cloud that we modeled. Assuming that the level ratio 
between the first fine-structure level and the Fe II ground 
state of this extra material is lower than 1/10, we estimate 
a lower limit to its distance of d = 9.5 kpc. 

Our successful fit of the UV pumping model to the ob- 
served evolution of the Fe II and Ni II excited levels has 
several interesting implications. 

The majority of the neutral gas closest to the GRB is 
at 1.7 kpc. If there would have been neutral material much 
closer in, we would not be able to reproduce the observed 
evolution of the excited levels with our model. Naturally 
this value is not accurate, mainly because of the uncertainty 
in the possible extinction in between the GRB and the 
cloud. However, this extinction is probably not very high. 
From the dust depletio n pattern in th e host (see Sect. [3] and 
IProchaska et al.ll2007l ; ISavaglioll200"6[ ). we estimate the ex- 
tinction to be low: Ay '^O.l. This value is in agreement 
with the host-galax y extinction estimate of Ay ~0.2 of 
in ison et ahl (|2006[ ). which is obtained from modeling the 
spectral slope. Moreover, the spectral slope from the UV 
pumping fit, P =— 0.5^° g is not very different from the ob- 
served spectral slope after correcting for the extinction in 
the foreground absorber at z = 1.1 {/3 = —0.8), while a 
large amount of dust extinction would severely affect the 
value of the spectral slope, provided that the extinction is 
not grey. 

We note that another lower limit to the absorb er dis- 
tance is set by the presence of Mg I (see [Prochask a et al.l 
l2006f ). assuming that it is at the same location as a 
large part of the Fe II and Ni II excited material. For 
GRB 051111. lProchaska et all (|2006l ) calculate that it Mg I 
would be at a distance smaller than 80 pc from the GRB, 
then Mg I woul d have been fully ionized. Using Eq. 2 of 
IProchaska et al.l (|2006t ). we estimate the lower limit to the 
distance where Mg I can survive to be c? = 45 pc for 
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GRB 060418. More recently. IChen et al.l ((2006) have esti- 
mated a lower limit to the distance from GRB 021004 to 
the absorbers along its sightline that are blue-shifted by 
2500-3000 km s"!. From the ratio of C II*/C II they find 
d > 1.7 kpc; limits for the other absorbers blue-shifted by 
less than 2000 km s^^ are not given. A similar distance 
limit is set for th e Si II gas associated with GRB 050730 
(jChen et al.ll2006[ ). All these distance limit estimates are 
fully consistent with our distance determination. 

The consequence is that any pre-GRB neutral cloud 
that was present at distances less than about 1.7 kpc, was 
severely affected by GRB 060418. Atomic species typical of 
the neutral ISM such as Fe II, Cr II, Zn II, etc., are likely 
ionized to a higher ionization stage. We note, however, that 
this does not imply that the GRB has ionized all neutral 
material up to this distance; it may be that the GRB ionized 
only its immediate surroundings, e.g. up to tens of parsecs, 
and that between the ionized region and 1.7 kpc, no signifi- 
cant amount of neutral material was present. But it is clear 
that the immediate environment of GRBs cannot be probed 
with these neutral ISM species, but possibly higher ioniza- 
tion lines may be detected. It is therefore of great interest 
to look out for higher ionization species not normally seen 
in optical spectra, that may originate from the immediate 
surroundings of the GRB. Possibly an ionization stratifica- 
tion could be observed, from higher ionization lines close 
to the GRB, to lower ionization species further out. X-ray 
spectroscopy instead of the optical will probably be the best 
tool to probe the immediate vicinity of the GRB. 

Because the photon energy threshold to ionize Fe II to 
Fe III is higher than the ionization potential of H I, this dis- 
tance limit also applies to neutral hydrogen, i.e. any signif- 
icant H I cloud that was present before the GRB exploded 
at distances smaller than approximately 1.7 kpc, will have 
been ionized. If we assume that GRB 060418 is not special 
with respect to other GRBs in this respect, most of the 
high-col umn density H I cloud s observed in GRB afte rglow 
spectra (jVreeswiik et al.|[2004l : iJakobsson et al.ll200(j ') may 
also be at typical kiloparsec distances. This assumption 
that GRB 060418 is not special, is supported by the lower 
limit to the distances of absorbers alo ng the GRBs 02100 4 
and 050730 sighthnes determined by IChen et al.l (|2006f) . 
which are also of kiloparsec scale. These H I clouds could 
either be part of a giant star-forming region in which the 
GRB was born, or simply clouds in the foreground in 
the host galaxy. We note that if the H I clouds are in- 
deed at kiloparsecs from the GRB, any metallicity esti- 
mate performed using optical spectroscopy is most likely 
not representative of the metallicity of the region where 
the GRB progenitor was born. The kiloparsec distance for 
these absorbers, combined with the significant differences 
between GRB-DLAs a n d QSO-DLAs in H I co lumn density 
( Vreeswi ik et al.' '2004| IJakobsson et al. 2006) , metallicity 
(jFvnbo et al... 2006 1 and dust depletion (Savaglio 2006. see 
also Sect. |3]), suggests that QSO sighthnes are not prob- 
ing the central kiloparsecs of (GRB) star-forming galaxies. 
This is consistent with our observation in Sect.jJl that the 
high-ionization profiles follow those of the low-ionization 
species very well in GRB 060418, while this is uncommon 
in QSO-DLAs. 

The gradual ionization of the neutral hydrogen close to 
the GRB could be observe d by monitoring the evolution 
of the Lya or metal line (jPerna fc Loebl Il998l we note 
that these authors suggest Mg II as metal line probe, but 



this line is normally highly saturated, especially in the high 
density GRB host galaxy environments, and therefore not 
suitable for this purpose). The redshift of GRB 060418 is 
too low for Lya to be covered by our UVES spectra. From 
the metal lines, we do not see any hint for such a gradual 
ionization; very likely observations have to be performed 
even quicker than the response time of our observations, to 
be able see this effect. 

Our UV pumping fit shows that it is possible to obtain 
the distance of the excited absorbing gas to the GRB. We 
have modeled the observations with only one cloud, but 
this can be extended to a multiple cloud model, where for 
each cloud one can obtain the distance with respect to the 
GRB, its velocity, its Fe II and Ni II column density and 
the cloud Doppler parameter, b. And for components not 
affected by the UV photons, a lower limit to the distance 
can be set. This way it could be possible to study the host- 
galaxy cloud structure, abundances and kinematics in more 
detail than before. 

Finally, the UV pumping fit not only constrains the 
properties of clouds in the GRB host galaxy, but also two 
properties of the GRB emission. One is the spectral slope 
of the GRB afterglow, even though this value is not very 
tightly constrained in our fit: /? =— O.Sj^^Q. The second 
property is the total UV flux that arrived at the cloud from 
the time of the burst trigger, i.e. optical flash (if present) 
and afterglow combined. This flux can be derived by de- 
termining the integral from fit parameter to to any time 
desired of the assumed light curve in the model. So even 
if no UV/optical observations were performed by robotic 
telescopes or Swift itself, the magnitude of a UV fiash 
can be constrained from a UV pumping fit. For the case 
of GRB 060418, we determine the limit on the total ob- 
server's frame V-band fiux that arrived from the GRB (UV 
flash and afterglow) from the time of the GRB trigger 
to the start of our flrst spectrum (i=11.4 minutes) to be 
(1.9 ± 0.3) X 10^^^ erg cm~^ Hz"^. For comparison, this 
flux is the same as that contained by a V=10 flash with a 
duration of 5 seconds. 

7. Conclusions 

Using the VLT Rapid-Response Mode in combination with 
UVES, we have obtained a unique time-series of high- 
resolution spectra of GRB 060418 with a signal-to-noise ra- 
tio of 10-20. These spectra show clear evidence for variabil- 
ity of transitions arising from the fine-structure levels of 
Fe II, and from metastable levels of both Fe II and Ni II. 
We model the time evolution of the Fe II and Ni II excited 
levels with three possible excitation mechanisms: collisions, 
excitation by IR photons only, and UV pumping. We find 
that the coUisional and IR photon scenarios can be rejected. 
Instead, the UV pumping model, in which a cloud with total 
column density N and broadening parameter 6 at a distance 
d from the GRB is irradiated by the afterglow photons, pro- 
vides an excellent description of the data. The best-fit val- 
ues are log iV(Fe U)^U.75toM, log ^(^i II)=13.84±0.02, 
d = 1.7 ± 0.2 kpc, and 6 = 25 ± 3 km s'^ The main 
consequence of this successful fit, is the absence of neu- 
tral gas, in the form of low-ionization metals or H I, at 
distances shorter than 1.7 kpc. Any pre-explosion neutral 
cloud closer to the GRB must have been ionized by the 
GRB. Therefore, the majority of very large H I column 
densities typically observed along GRB sighthnes may not 
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be in the immediate surroundings of the GRB; they could 
either be part of a large star-forming region, or foreground 
material in the GRB host galaxy. In either case, the metal- 
licity derived from absorption-line spectroscopy may not 
be representative of the metallicity of the region where the 
GRB progenitor was born. 
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